Unified Theory of Gliost Imaging with Gaussian-State Light 
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The theory of ghost imaging is developed in a Gaussian-state framework that both encompasses 
prior work — on thermal-state and biphoton-state imagers — and provides a complete understand- 
ing of the boundary between classical and quantum behavior in such systems. The core of this 
analysis is the expression derived for the photocurrent-correlation image obtained using a general 
Gaussian-state source. This image is expressed in terms of the phase-insensitive and phase-sensitive 
cross-correlations between the two detected fields, plus a background. Because any pair of cross- 
correlations is obtainable with classical Gaussian states, the image does not carry a quantum signa- 
ture per se. However, if the image characteristics of classical and nonclassical Gaussian-state sources 
with identical auto-correlation functions are compared, the nonclassical source provides resolution 
improvement in its near field and field-of-view improvement in its far field. 

PACS numbers: 42.30.Va, 42.50.Ar, 42.50.Dv 



I. INTRODUCTION 

Ghost imaging is the acquisition of an object's trans- 
verse transmittance pattern by means of photocurrent 
correlation measurements. In a generic ghost imaging 
experiment (see the example in Fig. [T]), a classical or 
quantum source that generates two paraxial optical fields 
is utilized. These fields propagate in two different direc- 
tions, through a linear system of optical elements that 
may include lenses and mirrors, and arrive at their re- 
spective detection planes. At one detection plane, the 
incident field illuminates a thin transmission mask, whose 
spatial transmissivity is the pattern to be measured, and 
is subsequently detected by a bucket detector that pro- 
vides no transverse spatial resolution. At the other de- 
tection plane, the incident field, which has never inter- 
acted with the transmission mask, is detected by a pin- 
hole detector centered at some transverse coordinate pi . 
The two photocurrents are then correlated and the out- 
put value is registered. This process is repeated as the 
pinhole detector is scanned along the transverse plane. 
The resulting correlation measurements, when viewed as 
a function of p^, reveals the power transmissivity of the 
mask. The image obtained by this procedure has been 
called a "ghost image," because the bucket detector that 
captures the optical field which illuminated the transmis- 
sion mask has no spatial resolution, and the the pinhole 
detector measures a field that never interacted with the 
transmission mask 

The first demonstration of ghost imaging utilized 
biphoton-state light obtained from spontaneous para- 
metric downconversion together with photon-counting 
bucket and pinhole detectors. This arrangement yielded 
a background-free image that was interpreted as a quan- 
tum phenomenon, owing to the entanglement of the 
source photons 0]. However, subsequent experimental 



d, 0] and theoretical [1, 0| considerations demonstrated 
that ghost imaging can be performed with thermalized 
laser light, utilizing either photon-counting detectors or 
CCD detector arrays to obtain ghost images, albeit with 
a background. 

The theory of biphoton ghost imaging requires quan- 
tum descriptions for both the optical source and its pho- 
todetection statistics, whereas thermal-light ghost imag- 
ing admits to a semiclassical description employing classi- 
cal fields and shot-noise limited detectors. This disparity 
has sparked interest [7j, l8|, |9|, [lOj in establishing a unify- 
ing theory that characterizes the fundamental physics of 
ghost imaging and delineates the boundary between clas- 
sical and quantum behavior. In this paper we develop 
that unifying theory within the framework of Gaussian- 
state (classical and nonclassical) sources. 

The foundation of our work is laid in Sec. [Til Here we 
begin by expressing the photocurrent cross-correlation — 
the ghost image — as a filtered fourth- moment of the field 
operators illuminating the detectors. Next, we briefly re- 
view the quantum-optics definition of classical states and 
specialize that discussion to zero-mean Gaussian states. 
Then, using the moment-factoring theorem for zero-mean 
Gaussian states, we obtain our fundamental expression 
for the Gaussian-state ghost image in terms of the phase- 
insensitive and phase-sensitive cross-correlations between 
the two detected fields, plus a background. The final 
part of Sec. [H] sets the stage for detailed understand- 
ing of ghost imaging by extending the standard theory 
of coherence propagation to include phase-sensitive field 
states. 
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Section lllll analyzes ghost imaging performed with 
three classes of Gaussian-state sources. We first con- 
sider a source possessing the maximum phase-insensitive 
cross-correlation — as constrained by its auto-correlation 
functions — but no phase-sensitive cross-correlation. Such 
a source always produces a classical state. Thermal light 
is of this class. We also consider a source with the maxi- 
mum classical phase-sensitive cross-correlation, given the 
same auto-correlations as in the previous case, but no 
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FIG. 1: (Color online) A simple ghost imaging setup. 



phase-insensitive cross-correlation. Finally, we treat the 
latter source when its phase-sensitive cross-correlation is 
the maximum permitted by quantum mechanics. The 
low-brightness, low-flux limit of this quantum source is 
the biphoton state. Thus these source classes span the 
experiments reported in 0, [3., JJ within a unified analyt- 
ical h'amework while admitting classical phase-sensitive 
light as a new possibility. In Sec. lIVI we discuss the image- 
contrast behavior that is obtained with these sources, and 
in Sec. El we generalize the Fig. [T] configuration to allow 
for a non-zero separation between the transmission mask 
and the bucket detector. We conclude, in Sec. IVII with 
a discussion of the ghost-imaging physics that has been 
revealed by our analysis. 



II. ANALYSIS 

Consider the ghost imaging configuration shown, using 
quantum field and quantum photodetection notation, in 
Fig. [TJ An optical source generates two optical fields, 
a signal Es{p,t)e~^'^"^ and a reference ER{p,t)e~^'^°\ 
that are scalar, positive frequency, paraxial field opera- 
tors normalized to have units photons /m^ s . Here, loq 
is their common center frequency and p is the transverse 
coordinate with respect to each one's optical axis. The 
commutation relations, within this paraxial approxima- 
tion, for the baseband field operators are [Tl| 



(1) 



[E„^{Pl,tl),El{p2,t2)] = d„^,iS{pi - P2)S{h - t2) , (2) 

where Sm,e is the Kronecker delta function, rrijt = S, i?, 
and d{-) is the unit impulse. Both beams undergo quasi- 
monochromatic paraxial diffraction along their respective 
optical axes, over an L-m-long free-space path, yielding 
detection-plane field operators 

Ee{p,t) = Jdp' Em{p',t-L/c)hL{p- p'), (3) 
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FIG. 2: (Color online) Photodetection model. 



where {£,m) = (1,5') or (2,i?), c is the speed of light, 
ko — LUo/c is the wave number associated with the cen- 
ter frequency, and h^ip) is the Huygens-Fresnel Green's 
function. 
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i2nL 



(4) 



At the detection planes, Ei(p,t) illuminates a quantum- 
limited pinhole photodetector of area Ai whose photo- 
sensitive region p e Ai is centered at the transverse 
coordinate pi, while E2{p,t), illuminates an amplitude- 
transmission mask T{p) located immediately in front of 
a quantum-limited bucket photodetector with photosen- 
sitive region p G A2- 

The photodetectors are assumed to have identical sub- 
unity quantum efficiencies and finite electrical band- 
widths, but no dark current or thermal noise (from sub- 
sequent electronics) contributes to the output current. 
Figure [2] shows the model utilized for the photodetec- 
tors, in which a beam splitter with field-transmissivity 
y/rj precedes an ideal photodetector to model the real 
detector's sub-unity quantum efficiency, and a low-pass 
filter with a real impulse response hg (t) follows the ideal 
photodetector to model the real detector's finite electrical 
bandwidth. It follows that the classical output currents 
from the two detectors correspond to the following quan- 
tum measurements [l^, [13, [I3| : 

im{t)^qjdu J dpE\ „^[p,u)E^^rn{p,u)hB{t - u), 

(5) 



for m — 1,2, where q is the electron charge. 



Er,,m{P^t) = VvEm{p,t) + \ - Evac,mip, t) , (6) 

and E^ac,m{Pi't) is a vacuum-state field operator. 

The ghost image at transverse location p^ is formed by 
time-averaging the product of the detector photocurrents 
to obtain an estimate of the ensemble-average equal-time 
photocurrent cross-correlation function, which is given by 

C{p,) = {h{t)h{t))^q\''A^ 

X / dp [dui [dU2hB{t~Ui)hB{t~U2)\T{p)\^ 

J A2 J J 

X {E\{p^,u^)El{p,U2)MPi,ui)E2{p.U2)) . (7) 
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where we have approximated the integral over the pin- 
hole detector's photosensitive region as the value of the 
integrand at times the photosensitive area, Ai. 

So far we have opted for the quantum description of our 
ghost imaging configuration, because it applies equally 
well to both classical-state and nonclassical-state sources. 
For the former case, however, we could have arrived at 
an equivalent answer by use of semiclassical theory. In 
particular, for a classical-state source we could replace 
the field operators with scalar classical electromagnetic 
fields, then employ scalar diffraction theory plus the shot- 
noise theory for photodetection to arrive at the pho- 
tocurrent correlation expression in Eq. ([7]), but with the 
field-operator fourth moment replaced by a classical-field 
fourth moment. Because a principal goal of this paper 
is to identify the classical/quantum boundary in ghost 
imaging, it is incumbent upon us to first review the quan- 
tum mechanical states that can generate photocurrent 
correlations which cannot be obtained from classical elec- 
tromagnetism and shot-noise theory. 



A. Semiclassical versus quantum photodetection 

Consider the ideal photodetector from Fig. [21 i.e, one 
with unity quantum efficiency, zero dark current and in- 
finite electrical bandwidth, for which individual photon 
detection events are registered instantaneously as cur- 
rent impulses carrying charge q. In semiclassical theory, 
the scalar optical field impinging on the photosensitive 
surface of the photodetector is a positive-frequency clas- 
sical electromagnetic wave, denoted by E{p,t)e~^'^°* . In 
accordance with the quantum description above, we as- 
sume that this field is paraxial, normalized to have units 
■\/ photons/m^s, and has center frequency ujq. Condi- 
tioned on knowledge of the field impinging on the pho- 
todetector, we have that i{t)/q, is an inhomogeneous 
Poisson impulse train with rate fmiction [l^, [l^ 

^i{t)^ f dp\E{p,t)\\ (8) 

J A 

where A is the detector's photosensitive region. Thus, 
regardless of whether the illuminating field is determin- 
istic or random, the photocurrent is subject to the noise 
that is inherent in this Poisson process, which yields the 
well known shot-noise floor of semiclassical photodetec- 
tion theory • Randomness in the illumination is then 
accounted for by taking E{p, t) to be a stochastic process, 
as is done in classical statistical optics [iTj . 

In the quantum theory of photodetection, the classical 
photocurrent produced by the same ideal photodetector 
is a stochastic process whose statistics coincide with those 
of the photon-flux operator measurement scaled by the 
electron charge [l^ . 

iit)^q [ ApE\p,t)E{p,t). (9) 

J A 



The photocurrent statistics are then governed by the 
state of the field operator t), so the shot-noise limit 
of semiclassical theory can be surpassed by some states, 
such as amplitude-squeezed states, or the eigenkets of 
continuous- time photodetection [ij, [3, [3 ■ 

In our quantum treatment of ghost imaging, the states 
of the optical field operator t) that we shall deem 
classical are those for which the measurement statis- 
tics predicted by quantum photodetection theory match 
those predicted by the semiclassical theory. It has long 
been known [l^, that when E{p, t) is in the coherent 
state |i?(p, t)), indexed by its eigenfunction E{p,t) and 
satisfying 

E{p,t)\E{p,t)) = E{p,t)\E{p,t)) , (10) 

the statistics of the i{t) measurement are identical to 
those from the semiclassical theory with the impinging 
classical field taken to be E{p,t). More generally, the 
two photodetection theories yield identical statistics for 
any quantum state that is a classical statistical mixture 
of coherent states — viz., for all states that have proper P- 
representations [l9| — when the classical field used in the 
semiclassical theory is comprised of the same statistical 
mixture of the coherent-state eigenfunctions [l^ IMI ■ 
Moreover, mixtures of coherent states are the only quan- 
tum states for which all quantum photodetection statis- 
tics coincide with the corresponding results found from 
the semiclassical theory. 

The quantum and semiclassical theories of photode- 
tection accord very different physical interpretations to 
their fundamental noise sources — quantum noise of the 
illuminating field versus shot noise arising from the dis- 
creteness of the electron charge — but for quantum states 
with proper P-representations their predictions are quan- 
titatively indistinguishable. So, because the correlation 
measurement at the heart of ghost imaging is a de- 
rived statistic from two photodetection measurements, 
the quantum theory of ghost imaging using states that 
have proper P-representations is equivalent to the clas- 
sical theory of ghost imaging using (classical) random 
optical fields plus shot-noise detection theory. Therefore, 
any truly quantum features of ghost imaging must be ex- 
clusive to optical field states that do not possess proper 
P-representations. 



B. Jointly Gaussian states 

Gaussian states offer both a practically relevant and 
a theoretically convenient framework for studying ghost 
imaging. Their practical relevance stems from thermal 
states and the biphoton state being special instances of 
Gaussian states. Their theoretical convenience arises 
from their being completely determined by their first 
and second moments, and from their closure under linear 
transformations. Moreover, as noted in Sec. I, Gaussian- 
state sources span the experiments reported in [2, [^, Q 
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and admit to the additional case of classical phase- 
sensitive light. Hence they provide an excellent unifying 
framework within which to probe the distinction between 
quantum and classical behavior in ghost imaging. 

Because the experiments in [2, S 01 employed 
zero-mean states, we shall assume that Es{p,t) and 
Eii{p,t) are in a zero-mean, jointly Gaussian state, 
i.e., the characteristic functional of their joint state 
has a Gaussian form [l^ specified by the (normally- 
ordered) phase-insensitive auto- and cross-correlations 
{El^{Pi,ti)Ei{p2,t2)), and the phase-sensitive auto- and 
cross-correlations {ErniPi,ti)Ee{p2,t2)) , where m^i ^ 
S,R. Because the experiments in [2|, H, 0] employed 
states whose phase-sensitive auto-correlations were zero, 
we shall assume that {E„i{Pi,ti)Em{p2i't2)) = for 
m — S, R. Finally, to simplify our analytical treatment, 
while preserving the essential physics of ghost imaging, 
we shall assume that the signal and reference fields are 
cross-spectrally pure, complex-stationary and have iden- 
tical auto-correlations, i.e., 

{El (p, , h )E^ {P2M))^ K^''^ (Pi , P2)^^"^ ih'h) (11) 
{^{p,,h)En{P2,t2))^Ki"Upi^p2)RsUt2-ti) (12) 
{Es{Pi,h)En{P2,t2))^Ki%iPi,p2)Ri%{h-ti) (13) 

for m — S, R, where the superscripts (n) and (p) label 
normally-ordered (phase-insensitive) and phase-sensitive 
terms, respectively. For convenience, and with no loss of 
generality, we shall assume that 

i?(")(0)=4'>) = 4^1,(0) = 1. (14) 

With the exception of the behavior of a background 
term, the physics of ghost imaging will be shown to arise 
entirely from the spatial terms in the preceding correla- 
tion functions. These will be taken to have Schell-model 
forms [Tgj . 

if(")(Pl,P2)=^*(PlM(P2)G(")(p2-Pl) (15) 

^S(Pi,P2) - A*ip,)AiP2)Gi%.iP2 - Pi) (16) 
4^kPi, P2) = A{p,)A{P2)Gi%,{P2 - Pi) , (17) 

with |A(p)| < 1, so that this function may be re- 
garded as a (possibly complex- valued) pupil function that 
truncates a statistically homogeneous random field with 
phase- insensitive auto-correlations G^"\p2 — Pi), phase- 
insensitive cross-correlation Gg"]j(p2 ^ Pi)i ^^'^ phase- 
sensitive cross-correlation Gg'^(p2 ~ Pi)- We shall also 
assume that G(")(p) is a real-valued even function of its 
argument :22]. Our task, in the rest of this subsection, 
is to establish the correlation-function bounds that dis- 
tinguish between classical and quantum behavior for the 
preceding jointly Gaussian states. 

Let us begin with Gaussian-state signal and refer- 
ence fields that have only phase-insensitive correlations. 



i.e., assume that {Es{Pi,ti)Eji{p2,t2)) = 0. Then, the 
phase-insensitive correlation spectra, given by the three- 
dimensional Fourier transforms 

5(")(fc,l]) = J^{g(")(p)^^"H^)} (18) 
~g^^},ikM) ^ .F{G^"],(p)4"],(r)} , (19) 

must satisfy [ij, [2^ the Cauchy-Schwarz inequality, 

\g^^Uk,n)\<~g^-\k,n), (20) 

from stochastic process theory |23|. Because the cor- 
relation spectra in (|20p fully determine the zero-mean, 
phase-insensitive, Gaussian state we are considering, 
this inequality is both necessary and sufficient to con- 
clude (via the equivalence developed in Sec. Ill Ap that 
all phase-insensitive Gaussian states have proper P- 
representations, and are therefore classical [2J|. The 
50/50 beam splitting of a continuous- wave laser beam 
that has first been transmitted through a rotating 
ground-glass diffuser — as was done in the experiments 
of d 13— yields signal and reference fields that are in 
a zero-mean, phase-insensitive, jointly Gaussian state in 
which (|20p is satisfied with equality. 

Now let us examine the more interesting case in which 
the zero-mean Gaussian-state signal and reference fields 
have a non-zero phase-sensitive cross-correlation, but no 
phase-insensitive cross-correlation. Here we will find that 
their joint state need not have a proper P-representation, 
viz., the state may be nonclassical. We have that the 
phase-sensitive cross-correlation spectrum of the signal 
and reference fields, 

g^^l,(fe,17) ^^{G%ip)Ri%ir)] , (21) 

satisfies [H 

\~9i%ik,n)\ < ^[l+5(")(fc,f])] g(»)(fe,r!), (22) 

whereas the Cauchy-Schwarz inequality for the phase- 
sensitive cross-correlation spectrum of a pair of classical 
stochastic processes imposes the more restrictive condi- 
tion [ll,[2fl] 

\g^^^^{k,n)\<~g(-Hk,n). (23) 

Zero-mean Gaussian states whose phase-sensitive cross- 
correlation spectra satisfy ([23|) have characteristic func- 
tionals consistent with that of a pair of classical stochas- 
tic processes. Hence these states have proper P- 
representations and are therefore classical. On the other 
hand, zero-mean Gaussian states whose phase-sensitive 
cross-correlation spectra violate have characteristic 
functionals that are inadmissible in stochastic process 
theory and are therefore nonclassical. In short, equal- 
ity in (|23p constitutes a well-defined boundary between 
classical and nonclassical zero-mean Gaussian states. 
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The difference between inequalities (|22)) and p3|) has 
a simple physical origin. Both derive from the fact that 
linear combinations of signal and reference fields have 
non-negative measurement variances. In the quantum 
case, however, the variance calculation leading to 122]) 
must invoke the field-operator commutators, whereas the 
derivation of P5)) has no such need. (Note that commu- 
tator issues do not arise in deriving (j20p . which is why 
this inequality is the same for the quantum and classical 
cases.) The upper bounds in ([^ and (^5]) are similar 
for 5^"^(fc,J7) 3> 1. Thus it might seem that there is lit- 
tle distinction between classical and quantum Gaussian 
states in this limit. While this will be seen below to be so 
for ghost imaging (when background is neglected), 50/50 
linear combinations of the signal and reference fields 
will be highly squeezed — thus highly nonclassical — when 
^("^(fe, n) > 1. At the other extreme, for g'-'^^k, fl) < 1, 
the quantum upper bound is approximately y/g^^^k, ft), 
which is significantly greater than the classical upper 
bound g^'^^k, fi). The phase-insensitive correlation spec- 
trum g^^'^k, fl) specifies the brightness of the signal and 
idler fields in units of photons. Thus g*^"^(fc, 51) ^ 1 is a 
low-brightness condition. In this regime we will see that 
there are appreciable differences between the ghost image 
formed with classical phase-sensitive light and quantum 
phase-sensitive light. 

Spontaneous parametric downconversion (SPDC), 
which was used in the original ghost imaging experi- 
ment [2], produces signal and reference fields that are 
in a zero-mean jointly Gaussian state with no phase- 
insensitive cross-correlation and no phase-sensitive auto- 
correlation, but with a phase-sensitive cross-correlation 
that saturates the upper bound in Furthermore, 
for continuous-wave SPDC operating at frequency de- 
generacy, this state is a two-field minimum-uncertainty- 
product pure state, generated by the Bogoliubov trans- 
formation [2^, [2^ 

Es{k,n) = AI{k,n)EsAk,^) 

+ vik,n)E],j-k,-n), (24) 

EB,{~k, -n) = AI{k, n)Eii^ (-fc, -n) 

+ Vik,n)El^{k,n) , (25) 

of the vacuum-state input fields, Es^ and En^ , where the 
transfer functions satisfy |M(fe,f7)p — |F(fc,ri)p = 1 to 
preserve the free-field commutator relations given in ^ 
and ©. 

In the low-brightness, low-flux regime, wherein 
5(")(fc,f7) < 1, \g<-P^k,n)\ « and at most 

one signal-reference photon pair is present in the electri- 
cal time constant of the ghost imager's photodetectors, 
the first and second-order moments of this Gaussian state 
match those obtained from the unnormalizable pure state 
comprised of a superposition of a dominant multimode 



vacuum and a weak biphoton component [25|, i.e., 

1^) - |o)s|o)fl + Jdk Jdn 

X (fc, n) e'^c^'^^) |fc, n)s\-k, ~n)R , (26) 

where 0(fe,fi) = Zg'^P'>{k,n) and |0)s|0)i^ is the muhi- 
mode vacuum state of the signal and reference. In the 
biphoton term, \k,n)s denotes the single-photon signal- 
field state with transverse wave vector fc and frequency 
detuning f2 from degeneracy; a similar interpretation ap- 
plies to the reference-field state | — fc, — fl)R. So, because 
the pure state given in Eq. ([26]) is the low-brightness, 
low-flux equivalent of the zero-mean jointly Gaussian 
state with maximum phase-sensitive cross-correlation 
and no phase-insensitive cross-correlation, it is clear that 
Gaussian-state analysis encompasses the previous bipho- 
ton treatments of ghost imaging using SPDC. 

As a final point about jointly Gaussian states, let us 
note how one may obtain classical phase-sensitive cross- 
correlations between the signal and reference fields. Such 
fields can be generated by imposing complex-conjugate 
zero-mean Gaussian-noise modulations, in space and 
time, on the fields obtained by 50/50 beam splitting 
of a continuous-wave laser beam. This saturates the 
upper bound in (|23p . because the resulting joint state 
is a Gaussian statistical mixture of the coherent states 
\E{p,t))s\E*{p,t))R. Existing modulator technology 
will limit the bandwidth achievable with such an ar- 
rangement to tens of GHz. Substantially broader band- 
widths might be realized by exploiting the classical (high- 
photon-flux) limit of nonlinear processes that generate 
phase-conjugate beams \2T\ . 



C. Coherence propagation 

The previous subsection laid out the statistical source 
models that we shall employ in our ghost imaging analy- 
sis; it was grounded in the second moments of the source- 
plane field operators Es and Er that completely charac- 
terize their zero-mean, jointly Gaussian state. However, 
our expression for the photocurrent correlation C{pi) in 
the Fig. [1] ghost-imaging configuration is given by Eq. (O, 
which requires a fourth moment of the detection-plane 
field operators Ei and E2- These detection-plane opera- 
tors result from L m free-space propagation of the source- 
plane operators, as given by Eq. ([3]). Jointly Gaussian 
states remain jointly Gaussian under linear transforma- 
tions, such as Eq. ([3]), and zero- mean states remain zero- 
mean as well. Thus, free-space diffraction over the L-m- 
long propagation paths transform the zero-mean, jointly 
Gaussian state of the source, with correlation functions 
given in Eqs. pTjl ^ P^ . into a zero-mean, jointly Gaus- 
sian state at the detection planes whose correlation func- 
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tions are cross-spectrally pure and given by 

{Ei{p„h)E,{p„t2))^K^^l{p„p,)R^f^l{t2~h). (28) 
In these expressions, 

<kPl, P2) - Jdp[ J dp', P2) 

^hl{p,-p[)hL{p2~p'2), (29) 
k[%„p,) = Jdp[ J dp', 4^],(pl, Pi) 

xhLip,-p',)hL{p2-p',). (30) 

for {m,m') = {1,S) or (2,i?), and likewise for {£,£'). 
Also, the temporal correlation behavior is unaffected by 
propagation, because the quasimonochromatic quantum 
Huygens-Fresnel principle, Eq. ([3]), only involves delay in 
time. It follows that the fundamental difference between 
the propagation of phase-insensitive and phase-sensitive 
correlation functions is the lack of conjugation in the 
propagation kernel of the latter, something which is re- 
sponsible for the propagation characteristics reviewed in 
Sec. Ed] 28]. 

Previous work on biphoton imaging has shown that 
the biphoton state propagates through free space in 
the same manner shown above for the phase-sensitive 
cross-correlation function [29| . This is not coincidental. 
We know from Eq. (pS)) that the biphoton wave func- 
tion is the phase-sensitive cross-correlation between sig- 
nal and reference fields with maximum phase-sensitive 
cross-correlation in the low-brightness, low-flux limit. 
For more general Gaussian states — which can have ar- 
bitrary brightness and photon flux and can be classical 
or nonclassical — it is necessary to consider the propaga- 
tion of the phase-sensitive cross-correlation function. 

Having related second moments of the detection-plane 
fields to their source-plane counterparts, we still need to 
find a fourth moment of those detection-plane fields in 
order to evaluate Eq. ([7]). For zero-mean jointly Gaus- 
sian states this ste p is easy. From the Gaussian moment- 
factoring theorem [l^ we find that the fourth-order mo- 
ment in Eq. ([7]) is given by 

{EliPi,ui)El (p, U2)Ei (pi, ui)i^2 (p, U2)) = 

{El (pi , Ui)Ei (pi , Ui))(El (p, U2)E2 (P, U2)) + 



where 



X i^i;)(Pi,Pi) / dp4';^{p,p)\T{p)\\ (33) 
is a non-negative non-image-bearing background, and 



C„ = q^Tj^Ai \r'[^> (t) P * he {t) * he {-t) 



\R^ll{t)\^*hB{t)*hBi~t) 



t=0 



(34) 
(35) 



are constants that depend on the temporal cross- 
correlations between Ei and E2- Here ★ denotes con- 
volution. 

The image-bearing term in C'(pi) is seen, from 
Eq. ((32)) . to be the object's intensity transmission pro- 
file, |T(p)p, filtered through a linear, space- varying filter 
whose point-spread function is given by a weighted sum 
of the squared magnitudes of the phase-insensitive and 
phase-sensitive cross-correlation functions at the detec- 
tion planes. In thermal-state ghost imaging, the phase- 
sensitive term vanishes, so that the point-spread function 
depends only on the phase-insensitive cross-correlation. 
In biphoton-state ghost imaging, the phase-insensitive 
cross-correlation is zero, thus yielding an image filter that 
depends only on the phase-sensitive cross-correlation. 
For general Gaussian-state signal and reference fields, 
however, both cross-correlations contribute to the image 
filter. 

Because the image-bearing part of Eq. p2|) only de- 
pends on the cross-correlations between the detected 
fields, whereas the non-image-bearing background de- 
pends only on the phase-insensitive auto-correlations, it 
is germane to note (see Appendix,) that any pair of phase- 
insensitive and phase-sensitive cross-correlation functions 
can be associated with a cZassicaZ zero- mean jointly Gaus- 
sian state, by appropriate choices of its phase-insensitive 
auto-correlation functions. Thus, if no constraint is 
placed on the background level in which the image is 
embedded, i.e., if the auto-correlation functions are not 
constrained, any image-bearing term attainable from 
Eq. ((32)) with a nonclassical Gaussian state source can be 
replicated identically by a classical Gaussian-state source. 
Hence, ghost-image formation is intrinsically classical. 



\{El{p,,U,)E2{p,U2))\' + m{Pl,Ul)E2{p,U2))\'m 

Substituting Eq. ((211) into Eq. ([7]), along with Eqs. dTTI) 
and (|28p , simplifies the photocurrent cross-correlation ex- 
pression to 

C(Pi) = Coip,)+Cn[ dp\K[';^ip„p)\'\Tip)\^ 

+ cj dp\K'f}{p,,p)\^\T{p)\\ (32) 

JA2 



D. Near- field versus far-field propagation 

Here we review the main results for paraxial, quasi- 
monochromatic, phase-insensitive and phase-sensitive co- 
herence propagation through free space [2^ that will be 
combined, in the next section, with Eq. ((32p to iden- 
tify the imaging properties of the Fig. [Tj configuration. 
Because Eqs. (I^Hl) and (150)) show that propagation only 
affects the correlation functions' spatial components, we 
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shall focus exclusively on them. In order to highlight the 
difference between the propagation of phase-insensitive 
and phase-sensitive coherence, we shall consider (real 
and even) Gaussian-Schell model spatial cross-correlation 
functions for both, i.e., we assume 30.] 



IP2-P1I /2po 



(36) 



for X — n,p, where P is the photon flux of the signal 
(and idler), ao is the attenuation radius of the trans- 
verse intensity profile, and po is the transverse coherence 
radius, which is assumed to satisfy the low-coherence con- 
dition Po ^ flo- 

We compare phase-insensitive and phase-sensitive cor- 
relation propagation in two limiting regimes: the near 
field, which corresponds to the region in which diffrac- 
tion effects are negligible, and the far field, in which 
diffraction spread is dominant. For phase-insensitive co- 
herence propagation, it is well known that a single Fresnel 
number, Dg = /coPoao/2i, distinguishes between these 
regimes, with Dq ^ 1 corresponding to the near field 
and Dq <^ 1 being the far field [3] ■ Note that this 
Fresnel number differs from that for the diffraction of 
a coherent laser beam with intensity radius ao, which is 
Dcoh = kQaQ/2L. This difference reflects the coupling be- 
tween coherence radius and intensity radius that occurs 
in free-space diffraction of partially-coherent light. In 
particular, far-field propagation of the phase-insensitive 
correlation function from Eq. (|36p results in an intensity 
radius satisfying ql — uq/Dq = 2L/A;oPo and a coher- 
ence radius given by pL = Po/Dq = 2L/kQao, i.e., the 
far-field intensity radius is inversely proportional to its 
source-plane coherence length and the far-field coherence 
length is inversely proportional to the source-plane inten- 
sity radius. 

The phase-sensitive correlation function from Eq. (|36p 
propagates in a distinctly different manner from its 
phase-insensitive counterpart. In this case we find that 
coherence-radius diffraction and intensity-radius diffrac- 
tion are decoupled 28] . Two Fresnel numbers are then re- 
quired to distinguish the near field from the far field: the 
Fresnel number for diffraction of the coherence length, 
Dn ~ fcoPo/2i; and the Fresnel number for diffraction 
of the intensity radius. Dp = fcoao/2L. The near-field 
regime for phase-sensitive correlation propagation occurs 
when both Fresnel numbers are much greater than one, 
and the far-field regime is when both are much less than 
one. Because we have imposed the low-coherence con- 
dition, Po ^ 0,0 J we can say that the near-field regime 
for phase-sensitive coherence propagation is Dn 3> 1 
and its far-field regime is Dp <C 1. Each of these con- 
ditions is more stringent than the corresponding condi- 
tion for phase-insensitive light. Nevertheless, the far-fleld 
propagation of the Gaussian-Schell model phase-sensitive 
correlation function from ([36]) still yields Po/Dq for the 
far- field intensity radius and ao/Do for the far-field co- 
herence radius. However, whereas the far-field phase- 
insensitive correlation is highest for two points with equal 



transverse-plane coordinates, the far-field phase-sensitive 
correlation is highest for two points that are symmetri- 
call y d isposed about the origin on the transverse plane 

Figure [3] highlights the difference between propagation 
of the phase-insensitive and phase-sensitive correlation 
functions from Eq. ((36)) . In this figure we have plot- 
ted the 

for X = 



e ^-attenuation isocontours for \K^^]^{pi, p2)\, 
n,p, versus the sum and difference coordinates 



Ps = (P2 + Pi)/2 and Pa 



P2 



All transverse- 



coordinate pairs that correspond to the interior region of 
a contour are both coherent and intense. From Eq. ([36]), 
it is straightforward to verify that all magnitude isocon- 
tours of our Gaussian-Schell model correlations are el- 
lipses. In the near field, because of our low-coherence 
assumption, the e~^-attenuation isocontours — for both 
the phase-insensitive and the phase-sensitive correlation 
functions — have their minor axes along the difference co- 
ordinate. In the far field, we find that diffraction leads 
to identical increases along the major and minor axes 
of the phase- insensitive correlation's e~^-attenuation iso- 
contour. For the corresponding far-fleld phase-sensitive 
correlation's isocontour we get inverted behavior, with its 
minor axis now aligned with the sum coordinate and its 
major axis along the difference coordinate. Thus, the far- 
fleld phase-insensitive correlation function is dominated 
by a narrow function in the difference coordinate \pa\, 
whereas the far-field phase-sensitive correlation function 
is a narrow function in the sum coordinate |p^|. 



4 |Pdl/2 



2L 



koPo 



2L 



koao 

po 



• Near field 

■ Far field, phase-sensitive 

■ Far field, phase-insensltlve 
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FIG. 3: (Color online) Isocontours corresponding to the 
e~^-attenuation levels for the phase-sensitive and phase- 
insensitive correlation functions in the near-field and the far- 
field regimes. 
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III. NEAR- AND FAR-FIELD GHOST IMAGING 
WITH GAUSSIAN-STATE LIGHT 

We are now fully equipped to compare the ghost- 
imaging performance achieved in the Fig. [1] configura- 
tion with various Gaussian-state sources. We shall as- 
sume that the signal and reference fields Es and Eji 
are in a zero-mean, jointly Gaussian state with identical 
phase-insensitive auto-correlations given by the following 
Gaussian-Schell model: 

iC(")(Pl,P2)i?(")(t2-tl) = 

where po <C ap, Tq is the coherence time, and P is 
the photon flux. We will begin our treatment with 
the thermal-state source, for which the signal and refer- 
ence have a non-zero phase-insensitive cross-correlation, 
but no phase-sensitive cross-correlation. As noted in 
Sec. IIIBl such states are always classical. 

A. Ghost imaging with phase-insensitive light 

Consider jointly Gaussian signal and reference fields 
with auto-correlations given by Eq. ([57]) and no phase- 
sensitive auto- or cross-correlations. Inequality ((^ 

implies that \{Eg{p-^^,ti)Eii{p2,t2))\ is maximum when 
it equals the auto-correlation function Eq. ([57)1 . We 
will take this to be so — to maximize the strength of 
the ghost image — and assume that this phase-insensitive 
cross-correlation is real-valued and non-negative. Be- 
cause near-field detection-plane correlations coincide 
with source-plane correlations, we can now obtain the 
near-field ghost image by substituting the right-hand side 
of Eq. ([57]) into Eq. (|32p . Doing so gives us the following 
result, 

C{p,) = Co(Pi)-f C„(2P/™g)2g-2lp,lVa^ 

X / dpe-^P^-P^^/P°e-^\P\^/''o\T{p)\^ . (38) 

Ja2 

Equation ((38|) reveals three significant features of the 
near-field, thermal-state ghost image. First, the ghost 
image is space-limited by the reference beam's average 
intensity profile, so that the object must be placed in the 
field of view ao jsH ]. Second, the useful transverse scan- 
ning range of the pinhole detector is restricted to the 
field of view oq. Finally, and most importantly, the finite 
cross-correlation coherence length po limits the resolu- 
tion of the image. When field-of-view limitations can be 
neglected, the ghost image in Eq. ([55]) is proportional 
to the convolution of the object's intensity transmis- 
sion, |T(p)p, with the Gaussian point-spread function 
g-lpl /Po. Thus the spatial resolution, defined here as 



the radius to the e ^-level in the point-spread function, 
is \/2/9o- 

Now let us suppose that the ghost image is formed in 
the far field, when Z^o <C 1, with the source correlations 
as assumed for the near-field regime. In this case we 
must first propagate source correlations — given by the 
right-hand side of Eq. (P7|) — to the detection planes via 
Eq. pS)) . It turns out that the detection-plane signal 
and reference fields still have maximum phase-insensitive 
cross-correlation, 

<l(Pl-P2)<^,(i2-il) = 

^e~(|pil'+|P2l')/4-|P2-PilV2p?,g-(t2-ti)V2To" (39) 

TTOI 

where m,£ G {lj2}, — 2L/kQpo and pL = 2L/fcoao, 
and the ghost image signature becomes, 

C(pi) = Co(pi)-)-C„(2PA4)2e-2|PilVai 

X / dpe-l''i-''l'/''^e-2|pl'K|r(p)|2_(4o) 

Ja2 

Therefore, the far-field field of view increases to 2L/kopo 
while the image resolution degrades to 2^/2L/koao, but 
the three conclusions drawn from the near-field image 
signature Eq. psp remain valid in the far-field regime. 
Because the resolution of the image degrades with prop- 
agation, so long as field of view is not the limiting factor, 
it is more desirable to place the object in the source's 
near field. 



B. Ghost imaging with phase-sensitive Hght 

Now we shall shift our focus to Gaussian-state signal 
and reference fields that have a non-zero phase-sensitive 
cross-correlation, but zero phase-insensitive cross- 
correlation. Applying the Cauchy-Schwarz bound ([^ to 
the Gaussian-Schell model auto-correlations in Eq. ((571) 
we find that the maximum \ {Es{pi,ti)Eii{p2,t2)) \ for a 
classical Gaussian state is also given by Eq. (|37p. Sim- 
ilar to what we did for the phase-insensitive case, we 
shall take the phase-sensitive cross-correlation to achieve 
its classical magnitude limit and assume that it is real- 
valued and non-negative. Then, because the detection- 
plane cross-correlation equals the source-plane cross- 
correlation in the near field, we can immediately get the 
near-field ghost image by substituting the right-hand side 
of Eq. daZl) into Eq. obtaining 

C{p,) = Co{p,) + Cp{2P/naire-'\f^\'/^o 

X / dpe-\P^-P\"/Poe-^\P\"/''o\T{p)\\ {41) 

JA2 

Equations ([M)) and ([55)1 give Cn = Cp, for our Gaussian- 
Schell model source, making the near-field ghost image 
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formed with classical phase-sensitive light identical to 
the near-field ghost image formed with phase-insensitive 
light. 

When the source-to-object separation is in the far- 
field regime for phase-sensitive coherence propagation, 
then the source-plane phase-sensitive cross-correlation 
that gave the preceding near-field ghost image gives rise 
to the following detection-plane phase-sensitive cross- 
correlation [28l. 



from which we see that the low-brightness regime cor- 
responds to PToPq/qI <^ 1. The source-plane phase- 
sensitive cross-correlation in this regime is then found to 
be 



{Es{Pl,tl)ER{p2,t2)) = (2/71 



,1/4 Oo 



^e~(l''il'+l''2l')/«Q-|P2-PilVpoe-(*2-*i)VTo" 



(46) 



2P 



irar 



g-dPil +IP2I )/<-|P2+Pil /2plg-(t2-ii)V2To'= 



which leads to 

C(pi)-C^o(Pi) + Cp(2PAai)2e-2|PilV4 



f dpe-|Pi+''l'/'''e-2|PlV-i|r(p)|2, 

J A, 



,(42) 



(43) 



for the far-field ghost image formed with classical phase- 
sensitive light. Again invoking Cp = Cn, for our 
Gaussian-Schell model source, we see that the far-field 
ghost image formed with classical phase-sensitive light 
is an inverted version of the corresponding far-field ghost 
image formed with phase-insensitive light, i.e., it has field 
of view and spatial resolution v^PL, as did the phase- 
insensitive ghost image, but the phase-sensitive ghost 
image is proportional to |T(— p)p ★ e"'''' Z''^- whereas 
the phase-insensitive ghost image was proportional to 

Finally, we turn to the ghost image produced using 
a nonclassical Gaussian state, i.e., one whose phase- 
sensitive cross-correlation violates In what follows 
we will restrict our attention to two limiting cases in 
which the phase-sensitive cross-correlation is coherence 
separable, so that we may utilize the machinery devel- 
oped earlier in this paper. In both cases we will take 
{Es{Pi,'ti)Eii{p2t2)) to be real-valued and non-negative 
with the maximum magnitude permitted by quantum 
theory. The limits of interest for this source will be those 
of high brightness, ^("^(fe,!]) ^ 1, and low brightness, 
5("^(fe,f2) ^ 1 when the source's auto-correlations are 
given by Eq. p7)) . 

At high brightness, the distinction between the cross- 
correlation functions of the quantum and classical phase- 
sensitive sources becomes insignificant, so that results 
given above for the ghost image formed with classical 
phase-sensitive light are excellent approximations for the 
quantum case. At low brightness, however, our assump- 
tions yield a phase-sensitive cross-correlation spectrum 
satisfying 



2(27r) 



1/4 



(44) 



(45) 



Note that (jiG]) is still a Gaussian-Schell, cross-spectrally 
pure correlation function, so that in the source's near 
field we get 



c(pi) = ao(pi) 



2 < ^ /2F 



-2|Pi|-/ao 



/ dpe-2|Pi-plVp;5e-2|plVa^|T(p)|2 

JA2 



(47) 



This near-field ghost image has the same field of view, 
ao, as the near-field ghost images formed with classical 
(phase-insensitive or phase-sensitive) light, but its spatial 
resolution, po, is a factor-of-\/2 better than the spatial 
resolutions of those classical near-field imagers. In ad- 
dition, the quantum case's image-to-background ratio is 
much higher than those of the classical imagers, because 
Qq/PToPq 1 in the low-brightness regime. 

The far-field ghost image for the nonclassical source 
is obtained by propagating its phase-sensitive cross- 
correlation from Eq. to the detector planes and sub- 
stituting that result into ([32]). The result we obtain is 



C(pi) = Co(Pi) 



P 



=-|PilV4 



dpe 



-Ipi+pI 



A2 



\np)? 



(48) 



Thus, the far-field resolution achieved with the quantum 
source equals those realized using the classical sources 
considered earlier, but the field of view has been increased 
by a factor of ^/2. It is worth pointing out that the 
quantum-enhancement factors — of spatial resolution in 
the near field and field of view in the far field — derive 
from the broadening of the weak spectrum, 5'^"'(fc,f2), 
when its square root is taken. That these enhancement 
factors both equal \/2 depends, therefore, on our choos- 
ing to use a Gaussian-Schell correlation model. Other 
correlation functions would lead to different enhancement 
factors. Finally, as found above for the near-field case, 
the quantum source yields dramatically higher image-to- 
background ratio in far-field ghost imaging than both 
its phase-insensitive and phase-sensitive classical coun- 
terparts. 
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IV. IMAGE CONTRAST 



with 



Thus far we have concentrated on the image-bearing 
terms in the photocurrent correlation from Eq. p2p . 
These image-bearing terms are embedded in a back- 
ground Co(pi), which, as we have seen in the preceding 
section, is much stronger for classical-source ghost imag- 
ing than it is for low-brightness quantum-source ghost 
imaging. It therefore behooves us to pay some atten- 
tion to the effect of background on ghost-imaging sys- 
tems. For the sake of brevity, we will limit our discus- 
sion to the near-field imagers; the far-field cases can be 
shown to have similar image-contrast issues. Also, we 
shall assume that the transmittance pattern being im- 
aged lies well within the field of view of all these ghost 
imagers, and restrict ourselves to considering the behav- 
ior of C{pi) in an observation region TZ that encompasses 
the image-bearing terms while satisfying \pi\ <^ oq. In 
this case 



C = 



maxTC[C(pi)] - miuTC [C(pi) 
Co(0) 



(49) 



is a meaningful contrast definition. Its numerator quan- 
tifies the dynamic range of the image-bearing terms in 
the photocurrent correlation C{pi), while its denomina- 
tor is the featureless background that is present within 
the observation region. 

For analytical convenience, let us take the baseband 
impulse response /is(i) to be a Gaussian with e^^- 
attenuation time duration T^, 



/^5(^) = e-«*^/^^^y^ 



(50) 



The contrast for the classical (phase-sensitive or phase- 
insensitive) ghost imagers then satisfies 

C(c) ^ Ci')C^t'\ (51) 
where the spatial (s) factor is given by 

^(c) ^ maXpjTe(Pi)] - minpjXc(pi)] 



L dp|T(p)|2 



(52) 



with 



IM^ I dpe-\P^-P\'''>l\T{p)\\ (53) 

being the point-spread degraded image of |T(p)p, and 
the temporal {t) factor obeys 



Ci^' = 1/v/l + (T,/2To)2 



(54) 



Likewise, for the low-brightness regime quantum imager 
we find that its contrast, C'-'-', factors into the product of 
a spatial term 

nil) - °o "^'^^Pi (Pi )] ~ minpi [^q (Pi )] f.r. 

' -''-^ Ldp\T{pW ' ^^^^ 



i|r(p)p, (56) 



being its point-spread degraded image of |T(p)p, times 
a temporal term 



Ci'^ = l/Jl + Ti/2TS. 



(57) 



The preceding classical and quantum contrast expres- 
sions possess interesting and physically significant behav- 
ior. We shall first explore the classical case. To get at 
its contrast behavior, we will assume that T{p) is a bi- 
nary amplitude mask, as has often been the case in ghost 
imaging experiments. It follows that 



where 



Ci^) « pI/At « 1, 
At^ Jdp\T{p)\^ 



(58) 



(59) 



and the inequality in ([55)1 holds because A-t/pq is ap- 
proximately the number of resolution cells in the ghost 
image. Combined with the fact that C^'^'' < 1, Eq. ((58|) 
shows that classical-source ghost imaging always has low 
contrast according to our contrast definition. This is 
why classical-source ghost imaging has been performed 
with thermalized laser light and has used AC-coupling 
of the photocurrents to the correlator [l^. Thermal- 
ized laser light is a narrowband source, for which Tgi <C 

(c) 

Tq so that Cf « 1. The use of AC-coupling implies 
that the correlator is estimating the cross- covariance be- 
tween the photocurrents produced by detectors 1 and 
2, rather than their cross-correlation. This ensemble- 
average cross-covariance is given by C{pi) — Co(pi), so it 
might seem that covariance estimation alleviates all con- 
cerns with the background term. Such is not the case. 
Even though the background term does not appear in 
the photocurrents' cross-covariance, its shot noise and ex- 
cess noise dictate that a much longer averaging time will 
be required to obtain an accurate estimate of this cross- 
covariance function, i.e., to get a high signal-to- noise ra- 
tio ghost image. Now suppose that classical-source ghost 
imaging is attempted using broadband light for which 
Td/To ~ 10"^, corresponding to a THz-bandwidth source 
and GHz electrical-bandwidth photodetectors. In com- 
parison with a narrowband classical-source ghost imager 
of the same photon flux P, the broadband imager must 
use a 10^-times longer time-averaging interval to achieve 
the same signal-to-shot-noise ratio. 

Turning now to the contrast behavior of the low- 
brightness quantum-source ghost imager, our assumption 
of a binary amplitude mask leads to 



PToAt 



> l/PTo 



(60) 
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thin- lens imaging theory, with the following results: 
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FIG. 4: (Color online) Ghost imaging setup with relay optics. 



because of our field-of-view assumption. Thus in broad- 
band, low-brightness, low-flux quantum ghost imaging 
we find that 



> i/PTrf > 1, 



(61) 



where the last inequality invokes the low-flux condition. 
This is why biphoton sources yield background-free ghost 
images 0,1111, despite SPDC being a broadband pro- 
cess. 



V. RELAY OPTICS 

Our analysis has assumed that the detector plane coin- 
cides with the object plane, but a realistic ghost-imaging 
scenario will likely require a separation between these 
two planes, as shown in Fig. HI In this flgure, the bucket 
detector is placed Lji m away from the object and we as- 
sume no control over this path, but we allow ourselves to 
freely modify the signal-arm path. Thus we place a focal- 
length-/ lens di m behind the object plane and d2 m in 
front of the detector plane, such that 1/di + l/d2 = 1//. 
In addition, because the optical path lengths may be dif- 
ferent, we introduce a (L/j — di — d2)/c post-detection 
electronic time delay to maximize the temporal cross- 
correlation of the two detected fields. The resulting pho- 
tocurrent cross-correlation is then 

C'(pi) = Co(Pi) + C„ / cip2|if(:t(Pi,P2)l' 

+ cj dp2\K'^P\,{p„p2)\\ (62) 

in terms of the phase-insensitive and phase-sensitive 
cross-correlations, k'^J^2'{Pit P2) fo'' — of the 

detected fields Ey and E2' ■ The magnitudes of these 
detection-plane cross-correlations are easily found from 



^S(Pl,P2) 



2ttLr 



K[':^\Mp„p')T{p') 



dp'e 



/-iko{2p2-p'-\p'f)/2Lf. 



(63) 



where M = ~d2/di is the signal-arm magnification fac- 
tor. For a sufficiently large bucket detector we can ap- 
proximate the integrals in (j62p as covering the entire 
plane, viz.. 



C„^ / dp2 \k[TI{Pi,P2) 



dP2\K[7:iip„P2)\' 



= M^Cr, 



dp\K["^,\Mp„p)\'\T{p)\\ 



(64) 



(65) 



for m = n,p, where we the last equality follows from 
Parseval's theorem. In this limit, C'{pi) — M'^C{M Pi), 
where C{pi) is given by ([5^ . Hence choosing di = 
d2 = 2/ will yield an inverted version of the object-plane 
ghost image. Image resolution and field of view are then 
determined by the phase-sensitive and phase-insensitive 
coherence properties of the object-plane fields, and the 
placement of the detectors relative to this plane only de- 
termines the signal-arm optics that are needed to obtain 
this object-plane ghost image. 



VI. DISCUSSION 

The fundamental source property that enables acqui- 
sition of a ghost image — whether the source is classi- 
cal or quantum — is the non-zero cross-covariance be- 
tween the photon-fiux densities of the two detected fields, 
i.e., the cross-correlation of the photon-fiux densities mi- 
nus the product of their mean values. In particular, 
the product of mean values generates the background 
term, while the cross-covariance produces the image- 
bearing terms. For zero-mean, Gaussian-state sources, 
the cross-correlation of the photon-fiux densities, which 
is a fourth-order field moment, reduces to a sum of terms 
involving second-order field moments. Consequently, 
both phase-sensitive and phase-insensitive field-operator 
cross-correlations can contribute to the ghost image. In 
the Appendix □ we will show that any pair of phase- 
sensitive and phase- insensitive cross-correlation functions 
can be obtained with two classical Gaussian-state fields, 
so long as there are no restrictions on these fields' auto- 
correlation functions. In this respect, the ghost image 
does not contain any quantum signature per se. However, 
if we compare sources that have identical auto-correlation 
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functions, we find that nonclassical fields with low bright- 
ness and maximum phase-sensitive cross correlation offer 
a spatial resolution advantage in the source's near field 
and a field-of-view expansion in its far field. The field 
of view in the near field and the resolution in the far 
field are determined by the beam sizes at the source, and 
hence are identical for classical and nonclassical fields. 

A number of recent publications have implied that 
ghost imaging with thermal-state light cannot be ex- 
plained by classical electromagnetic theory in combina- 
tion with semiclassical photodetection theory, but that 
a strictly quantum-mechanical interpretation involving 
nonlocality must be used to understand such experi- 
ments 10, |32j. A key conclusion from our paper, how- 
ever, is that the classical theory of ghost imaging is 
quantitatively indistinguishable from the quantum the- 
ory of ghost imaging for any optical source that is in 
a classical state, regardless of the propagation geome- 
try. Here, a classical state is one whose density opera- 
tor has a proper P-representation, so that its photode- 
tection statistics can be correctly quantified with classi- 
cal, stochastic-field electromagnetism and detector shot 
noise. Thermal light-whether it is broadband, such as 
natural illumination, or narrowband, such as thermal- 
ized laser light — falls precisely within this category of 
Gaussian states. Therefore experiments utilizing thermal 
light sources alone cannot validate the quantum descrip- 
tion. Furthermore, and perhaps more critically, there 
is no nonlocal interaction in thermal-light ghost imaging. 
In particular, because the joint state of the signal and ref- 
erence beams is classical — in the sense noted above — it 
cannot lead to a violation of the Clauser-Horne-Shimony- 
Holt (CHSH) inequality [11]. We reiterate that it is 
the non-zero cross-covariance between the photon-flux 
densities of the signal and reference fields that is re- 
sponsible for the image-bearing terms obtained from the 
Fig. [T] setup. For Gaussian-state sources, this detection- 
plane cross-covariance is found, by moment factoring, 
from the phase-insensitive and phase-sensitive field cross- 
correlation functions. These detection-plane field cor- 
relations follow, in turn, from propagation of the cor- 
responding source-plane field correlations through Lm 
of free space. Thus, two classical fields that are gener- 
ated in a correlated fashion at a source, yet propagat- 
ing paraxially in two different directions, will still ex- 
hibit spatio-temporal correlations on transverse planes 
that are equidistant from the source, even though these 
planes may be physically separated from each other. This 
concept is both well known in and central to classical 
statistical optics [l3, [13 • It is not at all related to non- 
locality in quantum mechanics, e.g., to violation of the 
CHSH inequality. 

It is worth connecting some of the analysis presented in 
this paper with recent theory for the coherence properties 
of biphoton wave functions, which has led to an elegant 
duality between the partial entanglement of biphotons 
and the classical partial coherence of phase-insensitive 
fields [295. As we have shown in Sec. IIIBl the bipho- 



ton state is the low-brightness, low- flux limit of the zero- 
mean jointly Gaussian state with zero phase-insensitive 
cross-correlation but maximum phase-sensitive cross- 
correlation. In this limit, the biphoton wave function 
is the phase-sensitive cross-correlation function between 
the signal and reference fields, and therefore the duality 
between phase-insensitive coherence propagation and the 
biphoton wave function propagation is rooted in the dual- 
ity between phase-insensitive and phase-sensitive coher- 
ence propagation [cf. Eqs. (|29p and ([50]) ]. Furthermore, 
classical fields may also have phase-sensitive coherence. 
Thus, to correctly understand the fundamental physics of 
quantum imaging, it is crucial to distinguish features that 
are due to the presence of this phase-sensitive correlation 
in the source fields from those that require this phase- 
sensitive correlation to be stronger than what is possible 
with classical (proper P-representation) states. The fol- 
lowing examples clearly illustrate our point. When ghost 
imaging is performed with phase-sensitive light, image 
inversion occurs in the far field for both classical and 
quantum sources. This inversion is entirely due to the 
difference between the free-space propagation of phase- 
sensitive and phase-insensitive correlations, and it is not 
necessary for the phase-sensitive coherence to be stronger 
than classical. On the other hand, the background-free 
nature of ghost images formed with SPDC light arises 
from that source's phase-sensitive cross-correlation be- 
ing much stronger the classical limit, as we showed in 
Sec.IVl 

In summary, we have used Gaussian-state analysis to 
establish a unified treatment of classical and quantum 
ghost imaging. Our analysis reveals that ghost-image 
formation is due to phase-sensitive and phase-insensitive 
cross-correlations between the signal and reference fields. 
Because arbitrary cross-correlations can be achieved by 
classical and quantum sources alike, image contrast is 
the only distinguishing feature between a source that is 
classical or quantum. In particular, we emphasize that 
a classical source with phase-sensitive cross-correlation 
can produce an identical image to that obtained with a 
biphoton source — up to a different contrast and hence 
signal-to-noise ratio — even for ghost-imaging configura- 
tions that utilize lenses, mirrors or other linear optical 
elements. If we compare ghost images from classical 
and quantum sources having identical auto-correlations, 
thereby fixing the background level, the low-brightness 
quantum source offers resolution enhancement in near- 
field operation and field-of-view enhancement in far-field 
operation, in addition to higher contrast in both regimes. 
Furthermore, because far-field spatial resolution and the 
near-field field of view are determined by source-plane 
beam size, they are identical for classical and quantum 
sources. Finally, the conclusions in this paper are not 
contingent on having coincident object and detection 
planes. They apply so long as the signal arm can be 
freely modified to transfer the object-plane correlations 
to the detection plane via an appropriately-positioned 
lens. 
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APPENDIX: CLASSICAL GAUSSIAN STATES 
WITH ARBITRARY CROSS CORRELATIONS 

Let us use Es{^) and Eb.{^) to denote the signal and 
reference field operators, where x — (p, t) conveniently 
combines their space and time arguments. In this ap- 
pendix we will construct a zero-mean, jointly Gaussian, 
classical state for these two quantum fields that has ar- 
bitrarily prescribed phase-insensitive and phase-sensitive 
cross-correlation functions, 



4p) 



i^^^jj(Xi,X2) = (i?5(xi)£;fl(x2)), 



(A.l) 
(A.2) 



respectively. We only require that both functions be suffi- 
ciently well behaved that they can be regarded as kernels 
which map the Hilbert space of square-integrable func- 
tions into itself. Under this regularity condition we can 
perform sin gula r-value decompositions of these continu- 
ous kernels [3J] to obtain 

oo 

i^^]j(xi,X2) = J2 VmC^*m(.^l)'i'm{^2) , (A.3) 
m— 1 

oo 

i^g']^(xi,X2) = /i^VmCxO^^Cxa) , (A.4) 

rn— 1 

where {^^(xi)}, {$„(xi)}, {■0m(xi)} and {^'^(xi)}, 
for 1 < m < oo, are four complete and orthonormal sets 
spanning square-integrable functions, and the coefficients 
r]„i and fim are real, finite, and non-negative for all m. 

Suppose we define two pairs of free-space, paraxial field 
operators, {i?5/(x),_Efl/(x)} and {Es"{'x.),Er"{x.)}, hav- 
ing the modal expansions 



£'s'(x) = as/^„j$,„(x) 

771 — 1 

OO 

Er'{x) ^ Y ai?',m0m(x) , 



m — 1 



and 



oo 

Ew'{x) = Y afl",m*m(x). 



(A.5) 
(A.6) 

(A.7) 
(A.8) 



In these expansions, {a^,m}, for i = S',S",R',R" and 
1 < m < oo, is a set of photon annihilation operators, 
with the canonical commutation relations lag. 



m 7 



^i,i'6m,m' and [a£,m, a^',m'] = 0. 

Now, let us put the modes associated with the {at^m\ 
into a zero-mean, jointly Gaussian state whose only non- 
zero phase-insensitive cross-correlations are. 



(A.9) 



and whose only non-zero phase-sensitive cross- 
correlations are 



whence 



(4,(X1)^;^,(X2)) =2if(;:]j(xi,X2) 

(i5"(xl)-Bi^"(x2)) - 2i^if]^(xi,X2) 



(A.IO) 

(A.ll) 
(A.12) 



Classical Gaussian states must have correlations that 
obey the Cauchy-Schwarz inequalities from ([^0]) and (|23p . 
see [20|. Thus Eqs. (jA.9p and (jA.lOp imply that the 
modal auto-correlations must obey 



2 

m 5 



and 



(A.13) 



(A.14) 



for 1 < TO < C)o. We will take the modal auto-correlations 
to equal these lower bounds, by assuming that 



''S".m 



as".r. 



''R/'m 



^R" .m} ^l^ra^ 



(A.15) 
(A.16) 



We shall also assume that all modal correlations — aside 
from those which have already been specified — vanish. 
Equations CO|) . COO)) . (|XT5)) and COB)) then de- 
termine the zero-mean, jointly Gaussian state of the 
four fields, (x), E'/j- (x), £^5" (x), E'^-/ (x)}, which is 
the tensor product of zero-mean, jointly Gaussian states 
of {£'s'(x),E_R'(x)} and {Es"{x),Er"{x)}, because all 
their cross-correlations are zero. Defining the signal and 
reference fields via 



Es{x) = {Es'{x)+Esn{x 
Er{x) = {Er,{x) + Er.{x 



(A.17) 
(A.18) 



thus yields a pair of field operators that are in a zero- 
mean, jointly Gaussian, classical state with the desired 
phase-sensitive and phase-insensitive cross-correlation 
functions. 
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